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AN EFFICIENT TREE DECOMPOSITION METHOD FOR 
PERMANENTS AND MIXED DISCRIMINANTS 

DIEGO CIFUENTES AND PABLO A. PARRILO 


Abstract. We present an efficient algorithm to compute permanents, mixed dis¬ 
criminants and hyperdeterminants of structured matrices and multidimensional ar¬ 
rays (tensors). We describe the sparsity structure of an array in terms of a graph, 
and we assume that its treewidth, denoted as w, is small. Our algorithm requires 
0{n2‘^) arithmetic operations to compute permanents, and 0{n^ + n3‘^) for mixed 
discriminants and hyperdeterminants. We finally show that mixed volume computa¬ 
tion continues to be hard under bounded treewidth assumptions. 


1. Introduction 
The permanent of a n x n matrix M is defined as 


Perm(M) : = En Mi, 


(d 


TT Z = 1 


where the sum is over all permutations tt of the numbers l,...,n. Computing the 


permanent is #P-hard 35 , which means that it is unlikely that it can be done efficiently 


for arbitrary matrices. As a consequence, research on this problem tends to fall into 
two categories: algorithms to approximate the permanent, and exact algorithms that 
assume some structure of the matrix. This paper lies in the second category. We 
further study related problems in structured higher dimensional arrays, such as mixed 
discriminants, hyperdeterminants and mixed volumes. 

The sparsity pattern of a matrix M can be seen as the bipartite adjacency matrix of 
some bipartite graph G. This bipartite graph fully encodes the structure of the matrix. 
We assume here that the treewidth a; of G is small. Even though it is hard to determine 
the treewidth of a graph, there are many good heuristics and approximations that justify 
our assumption 12 . We show an algorithm to compute Perm(M) in 0{n 2“) arithmetic 
operations. In this paper, the notation O ignores polynomial factors in uj. We note 
that the algorithm can be used over any commutative ring. 

The permanent of a matrix can be generalized in several ways. In particular, given a 
list of n matrices of size n x n, its mixed discriminant generalizes both the permanent 
and the determinant |4,21 . Our algorithm for the permanent extends in a natural way 
to compute the mixed discriminant. The natural structure to represent the sparsity 
pattern in this case is a tripartite (i.e., 3-colorable) graph. The running time of the 
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resulting algorithm is 0{'n?+n 3“^), where u is the treewidth of such graph. In particular, 
this algorithm can compute the determinant of a matrix in the same time. 

More generally, our methods extend to generalized determinants/permanents on ten¬ 
sors. A special case of interest is the multidimensional permanent [2,18,^. Another 
interesting case is the hrst Cayley hyperdeterminant, also known as Pascal determinant, 
which is the simplest generalization of the determinant to higher dimensions [5,14,25 


Note that unlike the determinant, the hyperdeterminant is ^^P-hard, in particular be¬ 
cause it contains mixed discriminants as a special case 21,22 . 

Given a set of n polytopes in M”, its mixed volume provides a geometric generalization 
of the permanent and the determinant 30 . We focus on the special case of mixed 
volume of n zonotopes. Although there is no “natural” graph to represent the structure 
of a set of zonotopes, we associate to it a bipartite graph that, when the mixed volume 
restricts to a permanent, corresponds to the matrix graph described above. This allows 
us to give a simple application for mixed volumes of zonotopes with few nonparallel 
edges. Nevertheless, we show that mixed volumes remain hard to compute in the general 
case, even if this bipartite graph has treewidth 1 and the zonotopes have only 3 nonzero 
coordinates. 

The diagram of Figure summarizes the scope of the paper. It presents the main 
problems we consider, illustrating the relationships among them. Concretely, an arrow 
from A to B indicates that S is a special instance of A. It also divides the prob¬ 
lems according to its difficulty, with and without bounded treewidth assumptions. In 
this paper we start from the simplest problems, i.e., permanents and determinants of 
matrices, moving upwards in the diagram. 



(this paper) 


Figure 1. Diagram describing the complexity relations of computing: 
determinants, permanents, mixed discriminants, hyperdeterminants and 
mixed volumes. 

The document is structured as follows. In Section we review the concept of 
treewidth and tree decompositions. We also present three graph abstractions of a 
sparse matrix. Among these graphs is the bipartite graph G described above, and a 
projection onto the column set. In Section we present a decomposition method. 
Algorithmic that computes the permanent based on the graph G^. We hrst use graph 
G^ because the decomposition algorithm is easier to explain in this case. In Section |C 
we extend this method to work with the bipartite graph G, as shown in Algorithmic 
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We provide a Matlab implementation of this algorithm. In Section we discuss the 
case of mixed discriminants, presenting a decomposition method for it. We also treat 
the case of generalized determinants/permanents on tensors. Finally, in Section]^ we 
discuss the case of mixed volumes of zonotopes. 

Related work. 


Permanents. The best known to date method for exactly computing the permanent of 
a general matrix was given by Ryser and its complexity is 0{n2^) [^. There are two 
main research trends on permanent computation: approximation algorithms and exact 
algorithms for structured matrices. We briefly discuss related work in both of them. 
Tree decomposition algorithms, which belong to the second group, will be presented 
afterwards. 

We hrst mention some work on approximation schemes. For arbitrary matrices, 
Gurvits gave a randomized polynomial time approximation, with error proportional to 
the operator norm 21 . For nonnegative matrices there is a vast literature, see e.g., 
and the references therein. Most remarkably, Jerrum et ah gave a fully polynomial 
randomized approximation scheme (FPRAS) [^. Recent work studies approximation 
schemes based on belief propagation, also for nonnegative matrices 


37,38 


As for exact algorithms for structured matrices, different types of structure have been 
explored in the literature. Fisher, Kasteleyn and Temperley gave a polynomial time 
algorithm for matrices whose associated bipartite graph is planar Barvinok 


showed that the permanent is tractable when the rank is bounded 
0/1 circulant matrices has also been considered 
matrices 


15 


The case of 
as well as sparse 0/1 Toeplitz 
Schwartz showed a 0(log(n)2^^) algorithm for certain band, Toeplitz 

Temme and Wocjan showed a 0(n2^"'^) 
Note that for arbitrary band matrices 


matrices, where w is the bandwidth 


31 


32 


algorithm for a special type of band matrices 
our algorithm is 0(n2^"'). 

Permanents and treewidth. Tree decomposition methods for permanent computation 
have been considered. Courcelle et al. hrst showed that the permanent can be com¬ 


puted efficiently if the treewidth is bounded 16 , although their methods, based on the 


Feferman-Vaught-Shelah Theorem, do not lead to an implementable algorithm. Later 


work of Flarup et al. gives a 0{n2^^‘^ algorithm 


20 


m 


This algorithm is extended 
26 to a wider class of matrices. Note the strong dependency on the treewidth. Fur¬ 


thermore, the graph abstraction used in the above methods, which is not the bipartite 
graph G, has two inconvenient features: its treewidth can be signihcantly larger than 


the one of G (see Example 2.1) and it is dependent on the specihc order of the columns 
of the matrix (see Remark 2.2). 

They gave a 0{n2‘^) 


Closer to this paper is the work of van Rooij et al. 


36 


decomposition algorithm for counting perfect matchings in a graph. Counting perfect 
matchings is closely related to the permanent, and one could derive from their proof an 
analogous, but different, method for calculating the permanent. Our algorithm could 
be seen as a variant of such method that is easier to extend to the higher dimensional 
problems we consider. 
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Mixed discriminants, mixed volumes, tensors. The higher dimensional problems we 
study generalize the permanent of a matrix, and thus are ^P-hard in general. As for 
the permanent, there are two natural relaxations: approximation algorithms and exact 
algorithms under special structure. Approximation algorithms have been considered 
for mixed discriminants [7,21 , mixed volumes and multidimensional permanents 

[^. As for exact algorithms under special structure, we are only aware of Gurvits’ 
tractability result for mixed discriminants and the 4-hyperdeterminant under some 
bounded rank assumptions 


21 


To our knowledge this is the hrst paper that studies tree decomposition methods 
for mixed discriminants, mixed volumes and generalized determinants/permanents on 
tensors. Related to this is a recent log-space algorithm for computing determinants 
under bounded treewidth assumptions [^. Also related is the problem of partitioning 


a low treewidth graph into fc-cliques, which is considered in 36 


2. Tree decompositions and matrix graphs 

In this section we review some basic facts regarding tree decompositions of a graph. 
We also present three graphs that can be associated to a sparse matrix. 


2.1. Tree decompositions. The notions of treewidth and tree decompositions are 
fundamental in many areas of computer science and applied mathematics 12,17 . Intu¬ 


itively, the treewidth of a graph G is a measure of how close it is to a tree. A graph has 
treewidth 1 if and only if it is a forest, i.e., a disjoint union of trees. The smaller the 
treewidth, the closer the graph is to a tree, and the easier it is to solve certain problems 
on it. We note that a graph of treewidth u has at most noj edges, and thus treewidth 
imposes a sparsity constraint. We give the formal dehnition now. 


Definition 2.1. Let G be a graph with vertex set X. A tree decomposition of G is a 
pair (T,x), where T is a rooted tree and y : T —)■ {0,1}^ assigns some y(f) C X to 
each node t of T, that satishes the following conditions. 

i. The union of {x(t)}t£T is the whole vertex set X. 

ii. For every edge {xi,Xj) of G, there exists some node t of T with Xi,Xj G y(t). 

iii. For every Xt E X the set {t : Xi E y(t)} forms a subtree of T. 

The sets y(t) are usually referred to as bags. The width uj of the decomposition is the 
size of the largest bag (minus one). The treewidth of G is the minimum width among 
all possible tree decompositions. 


Algorithms based on tree decompositions typically depend exponentially on the width 
of the decomposition, and polynomially on the number of nodes [^. Thus, given a 
graph G it is desirable to obtain a tree decomposition of minimum width. However, 
hnding the treewidth is NP-hard [^. The treewidth of some simple graphs are known: 
for a tree is 1, for a cycle graph is 2, for the n x n grid graph is n, for the complete 
graph Kn is n — 1, for the complete bipartite graph Kn,n is n. For general graphs, there 
are good heuristics and approximation algorithms 12 . 

Tree decompositions are closely related to chordal graphs 11,^. Indeed, given a 
chordal graph G, we can construct a tree decomposition where the bags y(t) correspond 
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to its maximal cliques. We remark that there are at most n maximal cliques in a chordal 
graph. In general, we can always assume that a tree decomposition has at most n nodes. 
The following is a simple property of tree decompositions. 

Lemma 1. Let (T, y) he a tree decomposition of G. Then for any clique Y of G there 
is some node t of T with Y C x(f). 

Proof. For each vertex y E Y, let Ty denote the subtree of all bags containing y. Let 
ty G Ty be the closest node to the root and let d{ty) denote the distance from ty to the 
root. Observe that if d{ty) < d{tyi) then ty' G Ty (otherwise, Ty P\Tyi = tj) and the edge 
(?/,?/') would not belong to any bag). Let t G {ty}y^Y be the farthest away from the 
root, i.e., d{ty) < d{t) for all y eY. It follows that Y C x(f). □ 


2.2. Matrix graphs. The sparsity structure of a matrix, i.e., its pattern of nonzero 
entries, can be described in terms of a graph. We consider here three possible graph 
abstractions of such sparsity structure, and we compare their treewidths. 

Let M he a n X n matrix. We will index the rows with a set A = {oi, ..., a^} and 
the columns with a set W = {xi, ... ,x„}. We use subindices to index the coordinates 
of M, i.e., Ma^x denotes the entry in the a-th row and x-th column. Similarly, let Ma be 
the a-th row of M. We now present two (undirected) graphs that are usually associated 
to a sparse matrix. 

Definition 2.2 (Bipartite graph). Let M be a n x n matrix, let A denote its set of 
rows, and let X denote its set of columns. The bipartite graph of M, denoted as G{M), 
has vertices AU X, and there is an edge (a, x) if Ma^x is nonzero. 

Definition 2.3. (Symmetrized graph) Let M be a n x n matrix, let A denote its set 
of rows, and let X denote its set of columns. The symmetrized graph of M, denoted as 
G^{M), has vertices 1,..., n and has an edge {i,j) if 7^ 0 or Ma^^xi 7^ 0. 

Remark 2.1. Note that G''^{M) is the adjacency graph of the symmetric matrix M+M'^ , 
assuming no terms cancel out. 

Remark 2.2 (Permutation invariance). Note that the permanent of a matrix is invariant 
under independent row and column permutations. The bipartite graph G preserves 
this invariance. On the other hand, the symmetrized graph G^ is only invariant under 
simultaneous row and column permutations. 


The bipartite graph G is our main object of study in this paper. Let u := tw(G) 
be the treewidth of G and Ug '■= tw(G^) the one of Tree decomposition methods 
based on graph G® have been studied before [3,16,20,26 . We claim that graph G 
is a better abstraction for the purpose of permanent computation. In particular, G 
preserves the permutation invariance of the permanent as stated above. Furthermore, 
Us can be much larger than u as shown in the following example. 


Example 2.1 (Two nonzero entries per row). Let M be a matrix with at most two 
nonzero entries per row. We claim that for all nontrivial cases the bipartite graph G 
has treewidth u <2. Let Gq be a connected component, and let hq be its number of 
row vertices. In order for Gq to have a perfect matching, it must have as many row 
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ai 

02 

03 

04 

O 5 

06 

07 

Os 

09 


Xi X2 X3 

-0 1 0 
0 0 1 
0 0 3 

0 0 0 
0 0 0 
0 3 0 

0 0 0 
0 0 0 
.3 0 0 


X4 X5 Xq 
2 0 0 
0 2 0 
0 0 2 
0 1 0 
0 0 1 
0 0 0 
0 0 0 
3 0 0 

0 0 0 


Xr Xs Xg 

0 0 0 - 

0 0 0 

0 0 0 

2 0 0 

0 2 0 

0 0 2 

3 1 0 

0 0 1 

0 0 0 . 



Figure 2. Graph abstractions of a matrix: bipartite graph G, sym¬ 
metrized graph G® and colnmn graph . 


vertices as colnmn vertices. Note also that Go has at most 2no edges becanse the row 
degrees are at most 2. Then Go is a connected graph with 2no vertices and at most 2 ?t,o 
edges, so it has at most one cycle. It follows that u < 2. 

On the other hand, we will see that Ug is nnbonnded. Let n = and consider the 
matrix M whose nonzero entries are 




M, 








M, 


^n — i 1 ^zm+1 


= 2 , 
. = 3, 
= 3, 


if m does not divide i 
ii 1 <i <n — m 
ii 1 <i <m 
ii 1 < i <m 


The corresponding graph G* contains the grid graph, and thns Ug > \/n. The case 
n = 9 is shown in Fignre 


The following example shows that the treewidth of G is always better than the 
treewidth of G*. 

Example 2.2 (G is “better” than G^). Let’s see that given a tree decomposition of G^ 
of width ojg we can form a tree decomposition of G of width 2ug. Let (T, i) be a tree 
decomposition of G®, where i{t) C {!,...,n}. Let /i : T —)■ {0,1}'^^'’^, be snch that 
/r(t) = {tti : i G i(t)} U{a;j : i G Then (T,/i) is a decomposition of G of width 2ujg. 

On the contrary, for a hxed u the treewidth of G^ is nnbonnded as seen in Example |2.1[ 

We now introdnce a third graph that we can associate to matrix M. 

Definition 2.4. (Colnmn graph) Let M be a n x n matrix, let A denote its set of rows 
and let X denote its set of colnmns. For any a G A let X (a) denote the set of nonzero 
components of row Mq. The column graph G^{M) has X as its vertex set, and for each 
a E A we pnt a cliqne in X{a). Eqnivalently, there is an edge {xi,Xj) if there is some 
a E A snch that Xi,Xj E X{a). 

Graph G^ can be seen as a projection of G onto the colnmn set X. We show in the 
following example that u < ojx + 1, where ux '■= tw(G^). 
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Example 2.3 {G is “better” than G^). Let {T,x) be a tree decomposition of G^ of 
width ux- For each row a & A we associate to it a nniqne node t E T snch that 
X{a) C x(f)- This assignment can be made becanse of Lemma For some t E T, 
let a{, be all rows that are assigned to t. Let’s replace node t of T with a 

path ti,t2 ■ ■ ■, tk, and let = x(^) U {a*} for j = 1,..., k. The nodes previonsly 

connected to t can be linked to any of the new nodes. By doing this for every t E T, 
we obtain a tree decomposition (T, p) of G of width ux + 1- 

On the other hand, for a hxed u the treewidth of G^ is nnbonnded. For instance, 
consider the matrix M whose only nonzero entries are: Ma-^xi = 1 (diagonal) and 
Mai,xi = 1 (first row) for all i. In this case G^ is the complete graph {ujx = n — 1) and 
G is a tree {u = 1). 

The reason why we consider the graph G^ is that we can give a very simple algorithm 
for the permanent based on it. We present this algorithm in Section]^ and then extend 
it to G in Section HI 

To conclnde this section, let’s see how the three graphs G, G"^, G^ can captnre the 
special case of a band matrix. 


Example 2.4 (Band matrix). Let Wi,W 2 E Z>o and let M be snch that Ma^^xj = 0 if 
either i — j > wi or j — i > W 2 . Let T be a path ti, ^ 2 ,0,.... Optimal decompositions 
(T, x), (T, l), (T, p) for the graphs G^, G^, G respectively, are given by 


X{ti) := {xi,Xi+i,... i{ti) := {z, i + 1,..., 7 + max{M;i, Wa}}, 

/^(Oi—l) • {(^i+wiy Xi, Xi - 1-1 , ... , Xj-|-jux +«;2 — 1 }) ■ '{ai+ioi 5 Xi-^.l, Xi-^-l , . . . , } 


The widths of these decompositions are u = uJx = wi + W 2 and ujs = maxjtci, W 2 }- 


3. Column decompositions 


Notation. For sets Y,Yi,Y 2 , let Y = Yi U 1^2 denote a set partition, i.e., Y = Yi U 1^2 
and hi n >2 = 0- 


In this section, we develop an algorithm to compnte the permanent based on a tree 
decomposition of the colnmn graph G^ (see Definition 2.4). For this section only, we 
denote the treewidth of G^ by uj . We will show that we can compnte Perm(M) in 
0(n4“’). Recall that the notation O ignores polynomial factors in uj. 

As before, A denotes the row set and X the colnmn set. We nse snbindices to index 
the coordinates of M, i.e., Ma^x denotes the element in row a and colnmn x. The 
following example illnstrates the methodology we follow. 



Figure 3. Matrix M, its colnmn graph G^ and a tree decomposition T. 
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Example 3.1. Consider the 5x5 matrix M of Figure [^and the following partition of 
its rows: 


A = {ai, 02 } U { 03 } U { 04 , as}. 


There is a simple expansion of Perm(M) in terms of this partition: 


Perm(M) = Perm 


Perm 


Perm 


Mat 

,Xi 

Mat 

Ma, 

,Xi 

Ma, 

'Mat 

XI 

Mat 

Ma, 

XI 

Ma, 

'Mat 

XI 

Mat 

Ma, 

XI 

Ma, 


Perm ( [Ma^^x ^]) Perm 
Perm ([Mag^^ 3 ]) Perm 
Perm () Perm 


M, 

M 


a4,a?2 


a5,X2 


M, 

M, 


a4,X2 


as, 2:2 


M, 

M, 


a 4 ,X 3 


as, 2:3 


M, 

M 


a4,tC5 


a5,tC5 


M, 

M, 


a4,X5 


as,a:s 


M, 

M, 


a4,X5 


as,a;s 


This expansion implies that to compute Perm(M) we just need to evaluate two 2x2 
permanents corresponding to {04,02}, three 1 x 1 permanents corresponding to {03}, 
and two 2x2 permanents corresponding to (04, 05 }. This requires only 14 multiplica¬ 
tions, compared to 4 x 5! = 480 multiplications using the dehnition. The reason why 
this formula exists is because the column graph of matrix M has a simple tree 
decomposition, which is shown in Figure 


As in the example above, we can always obtain an expansion of Perm(M) using a 
tree decomposition of graph . By carefully evaluating this formula we will obtain a 
dynamic programming method to compute Perm(M). 


3.1. Partial permanent. In our notation, the permanent of M can be expressed as 

Perm(M) = En 

TT asA 

where the sum is over all bijections vr : A —)■ X. For a given row o, let X(a) denote the 
column coordinates where it is nonzero. We will refer to a bijection tt as a matching if 
7 r(a) G X(a), i.e., Ma^n{a) 7 ^ 0, for all a & A. Then we can restrict the above sum to be 
over all matchings. 

We consider a tree decomposition (T, y) of the column graph G^. Note that by 
construction of G^ then X(a) is a clique for any a E A. Thus, Lemma says that we 
can assign each row o G A to some node t, such that X(a) C y(t). From now, we fix 
an assignment of each o to a unique node. Let At denote the rows that are assigned to 
node t. Thus, we have the partition A = UteT^i- 

For a node t G T, we denote the subtree of T rooted in t by Tt. Let y(Tt) be the 
union of x{t') over all t' E Tt, and similarly let At^ be the union of Af over all t' E Tt. 
For instance, for the root node we have Ay^^j = A and xiTroot) = X. 

For a hxed matrix M and some sets D T A and X C X we denote 

(1) perm{D, Y) := En 

TT a^D 

where the sum is over all matchings n : D ^ Y. Equivalently, it is the permanent 
of the submatrix of M corresponding to such rows and columns. Clearly, this only 
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makes sense if \D\ = |y|, and otherwise we can define perm{D,Y) = 0. We refer to 
perm{D,Y) as a partial permanent. 

3.2. Decomposition algorithm for the permanent. Algorithm [T] presents onr dy¬ 
namic programming method to compnte Perm(M). We will explain and derive this 
algorithm in the following sections. For each node t of the tree, the algorithm compntes 
a table Pt indexed by snbsets Y of x(t). It starts from the leaves of the tree and recnr- 
sively compntes the tables of all nodes following a topological ordering. The permanent 
of M is fonnd in the table corresponding to the root. 


Algorithm 1 Permanent with colnmn decomposition 

Input: Matrix M and tree decomposition (T, x) of colnmn graph G^{M) 
Output: Permanent of M 
1: procedure ColsPerm(M, T, x) 

2 : assign each a G A to some t with X{a) C x{t) 

3: order := topological ordering of T starting from its leaves 

4: for t in order do 

5: Qt := SUBPERMS(t, M) 

6 : if t is a leaf then 

7: Pt := Qt 

8 : else 

9: Cl,..., Cfc := children of t 

10: Pt :=EVALRECURSION(f, Qt, Pc^,..., PcQ 

11: return Proot{x{.root)) 

12: procedure SuBPERMs(f, M) 

13: At := rows assigned to node t 

14: Qt{y) '■= perm{At,Y) for all Y C x{t) 

15: procedure EVALRECURSION(t, Qt, Pci, • • •, PcQ 
16: for Cj child of t do 

17: Aj := x(cj) \ x(f), Aj := x{cj) f] xiQ 

18: Qc, (Y) := Pc^. (F U A,) for all F C Aj 

19: Pt := SUBSETCONVOLUTION(Qi,Qci,. • -yQcQ 

20: procedure SubsetConvolution(Po, Pi, ..., PQ 
21: P{Y):= Po{Yo)Pi{Yt)---Pk{YQ 

you-uYfc=y 


Algorithm has two main rontines: 

• For any node t, SubPerms compntes a table Qt with the permanents of all 
snbmatrices corresponding to t. (See Section 3.3). 

• EvalRecursion compntes table P* * of an internal node t, by combining table 
Qt with the tables Pc^,... ,Pcf. of the node’s children. (See Section 3.4). 
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The values Pt{Y) that we compute correspond to a partial permanent of the matrix, 
as we explain now. Consider the collection 

5:= {Y:x{Tt)\x{t)YYCx{T,)}. 

Observe that F G S' is completely determined by F fl xi^)- Therefore, if we let Y := 

Y nx(t), there is a one to one correspondence between S, and the collection S := {Y : 

Y C y(t)}. Then the partial permanents that we are interested in are 

Pt{Y) := perm{ATt, Y) = perm{AT„Y U {x{Tt) \ x{t)))- 


The reason why we index table Pt with Y instead of F, is that in this way it becomes 
clearer that the recursion formula is actually a subset convolution. 

Observe that the permanent of M is indeed computed in Algorithm [T| as for the root 
node we have Prootixii^oot)) = perm{A, X) = Perm(M). Also note that for a leaf node 
we have Pt{Y) = perm{At,Y) = Qt{Y). 


Example 3 . 2 . Consider the matrix M and tree decomposition T of Figure Let 
ti,t 2 it 3 be the nodes of T, where the central node ^2 is fhe root. We show the tables 
computed by Algorithm [T] The tables Qt with the permanents of all submatrices are: 


QH{{x,y}) = perm{{ai,a2},{x,y}), 
Qtiiix^y}) = perm{{a4,a5},{x,y}), 
Qt2i{x}) = perm{{a3}, {x}). 


for x,y e xih) = {xi, X3, x^} 
for x,y e xih) = {2^2, X 3 , X 5 } 
for X e xih) = {x2,X3,X4} 


We now show the final tables Pt for each node. For the leaves ti, fs we have Pt^ 
Pts = Qts- As for the root ^ 2 , the recursion is: 


Qti ) 


Pt2{{x2,X3,X4}) = Qt^{{x4})PtA{Xl,X3})Pt3{{x2,X5}) + 

Qt2{{x3})Ph{{xi, X4})Pt3{{x2, X5}) + Qt2{{x2})PtA{xi, X4})Pt3{{x3, X5}). 
Note that this recursion matches the permanent expansion in Example |3.1[ 

In the following sections we explain the two main routines of Algorithm [T| i.e., Sub- 
Perms and EvalRecursion, obtaining complexity bounds for them. 


3.3. Permanent of all snbmatrices. Let Mq be a rectangular matrix with row set Aq 
and column set Xq. As a part of our algorithm, which can be seen as the base case, we 
require a good method to compute the permanents of all submatrices of Mq. In other 
words, we want to obtain the partial permanents perm{D,Y) for all pairs {D,Y). We 
can do this in a very simple way using an expansion by minors. The following lemma 
explains such procedure and gives its running time. 

Lemma 2 . Let Mq be a matrix of dimensions ni x 77,2. Let Aq denote its row set, Xq its 
column set and let S = {{D, Y) Aq x Xq : \D\ = |E|}. We can compute perm{D, Y) 
for all {D,Y) E S in where Umax = max{ni,n2}. 

Proof. Let St = {{D, Y) : |T*| = |E| = i] for 1 < i < min{ni, 712 }. We use an expansion 
by minors to compute perm{D, Y) for {D, Y) G St, using the values of St^i. Let Qq be 



11 


the first element in D, then 


perm{D, Y) = V M. 


x^Y 




perm{D \ao,Y \x). 


Thus, for each [D, Y), we loop over at most n 2 elements, and for each we need 0{n. 
to find the sets -D \ Oq and Y \x. The result follows. 


max ) 
□ 


3.4. Recursion formula. The heart of Algorithm [T] is given by the recursion formula 
used, i.e., the procedure to obtain table P* of node t from the tables of its children. 
This recursion formula is given in the following lemma. 

Lemma 3. Let M be a matrix with associated column graph . Let (T, x) a tree 

decomposition of . Let t be an internal node ofT, and let Y be such that 

(2) xm)\x(t)CFCy(T*), \Y\ = \At,\ 

Let Cl,..., Cfc be the children of t. Then 

k 

(3) perm{ATt, ^ perm{At, Yt) perm{AT ,., Y^.) 

y i=i 

where perm{-, ■) is as in ([^ and the sum is over all y = (Yt, Fci, • • •, Yc^) such that: 

(4a) Y = YtU{Y,,U---UYj 

(4b) x(T,J \ x{t) C C xiY,) Yt C x{t). 

Proof. Observe that can be partitioned as 

Axt = AtU {At^^ U ■ ■ ■ U At^J. 

Let TT : —?• F be a matching. Let c be a child of t and let tTc : At^ —)■ F be the 

restriction of tt to At^. Let Yt := vr(Ai) C x{t) be the range of tt restricted to At, and 
let Fc be the range of tTc. As tt is injective, then equation (4a) holds. Observe also that 


Fc = 7r(ATj F x(Tc). Note now that if x G x{Tc) \ x{t)i then it is in the range of tt. 
However, as x ^ x(f) then x ^ x(T'c') for any other child c', and thus x has to be in the 
range of tt^. Thus, the range of tTc, i.e., Fc, contains x(^c) \ xif)- 

Therefore, for any matching tt : A^^ —)■ F and for any child c, vr induces a matching 
from At^ to some Fc that satisfy equations (|^. On the other hand, given Yt,Yc^,... 
satisfying (|^ and matchings tt^, tTc^, ... on At, At^^ ,... with such ranges, we can merge 
them into a function on Aj-^. Observe that (|^ ensures that the ranges of these matchings 
are disjoint and their union is F. We conclude that 


perm{A'rt,Y) = E n Ma, 


7r(a) 




E E 


at,nt(at) | ( J_ J_ (a^) 

Cj V flc 


y TTtiAt^Yt \ at 

TTc-Axq—^Yc 


^ perm{At, Yt) perm {At,. , Yc.). 

y c. 
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□ 


At first sight, the recursion of equation ([^ looks difficult to evaluate. It turns out 
that this formula is a subset convolution and thus it can be computed efficiently using 


the algorithm from 10 , as explained in the following lemma. We follow this approach 


in method EvalRecursion of Algorithm 
Lemma 4. Given the values of the partial permanents perm{At, Yt) and perm{AT^ , Yc^), 
we can evaluate equation ([^ for all Y satisfying ([^ in 0{k2^). 

Proof. Let Y := Y D x{t), Yt := Yt, Y^.. := W, H x(t), and let 

At := x(Tt) \ xif) Ac, := x{Tc^) \ xif) 

Pt(Y) := perm{ATt, E U At), QtiXt) := perm{At, Yt), Qc, (W,) := perm{AT,., Y^. U A^J 
Then equation can be rewritten as 


(5) 


P.{Y)= E QAn)J{Qc,(X,) 

rtuRiU-uRj^=y j=i 


where Yt Y x{t) and Yc. C xi^j) H x(t). Equation (|^ is a subset convolution over the 
subsets of x(t)- Therefore, it can be evaluated in Oikw"^ 2'“), where w = |x(t)|, using 
the algorithm from 


10 


□ 


The following theorem gives the running time of Algorithm proving that we can 
efficiently compute the permanent given a tree decomposition of of small width. 


Theorem 5. Let M be a matrix with associated column graph . Let {T,x) be a tree 
decomposition of of width oj. Then we can compute Perm(M) in 0{n4‘^). 

Proof. Let t be some node in T. We compute perm{ATt,Y) for every Y satisfying (|^. 
In particular, we will obtain Perm(A) = perm{AT,.„^t,x{Troot))- We will show that for 
each t we can compute perm^Axt, Y) for all Y in 0{{kt + 1)4^^), where kt is the number 
of children of t. Note that Ylt number of nodes of tree T. As the tree has 0{n) 

nodes, the total cost is then 0{n4P), as wanted. 

The base case is when t is a leaf of T, so that = A* and x(Tt) = x(f)- Let Mq 
be the submatrix of M with rows At and columns x(t). Then all we have to do is to 
obtain the permanent of some submatrices of Mq. Observe that \At\ < |x(^)| A as 
otherwise there is no Y satisfying (§. Thus, we can do this in 0(2^^^) using Lemma |2[ 

Assume now that t is an internal node of T with kt children and let Y that satishes (^. 
Then equation tells us how to hnd perm{ATt,Y). Lemmaj^says that we can evaluate 
the formula in 0{kt 2‘^), assuming we know the values of the terms in the recursion. Note 
that we already found perm^Ax^, Yc) for all children and that we can hnd perm{At, Yt) 
for all possible T) in 0(2^^) in the same way as for the base case. Thus, it takes 
0(4^ + kt 2^) = 0{{kt + 1)4^^) to compute perm{ATt, Y) for all Y. □ 


Remark 3.1. The factor 0(4^^) in the proof came from the base case. This bound can 
be improved, but we omit this as for the approach of Section]^ based on the bipartite 
graph, the bound will be 0{2‘^). 
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3.5. Computing the determinant. Given the similarity between permanent and 
determinant, it should be possible to hnd an analogous algorithm for the determinant. 
We will derive such algorithm in this section. Ironically, this algorithm is slower than 
the one for the permanent. The reason is that the approach we follow does not take 
advantage of linear algebra: we loop over all permutations (carefully) and then compute 
its sign. We remark that our algorithm does not use divisions and thus can be applied 
in any commutative ring. The ideas from this section will be used in Section to derive 
a decomposition algorithm for the mixed discriminant. 


Example 3.3. Consider again the matrix M of Figure and observe that a similar 
expansion holds for the determinant: 


Det(M) = Det 


— Det 


Det 


M, 


aiyXi 


M, 

M, 


ai,xi 


a2,xi 


M, 

M 


ai,xi 


a2,xi 


M, 

M, 


ai,xz 


^2,^3 


Ma 


M, 

M 


ai ,X 4 




Det Det 

Det {[Ma 3 ,xi]) Det 
Det ([Mag , 0 , 2 ]) Det 


M, 

M 


0.4 ,^2 


CLb^X2 


M, 

M, 


M, 

M, 


M, 

M 


a4,a:5 


^5,3^5 


a4,a;2 


as ,^2 


M, 

M 


a4 


as,2:5 


a4,X3 


a5,X3 


M, 

M 


a4 ,xs 


as 5 X 5 


As suggested in the above formula, the recursion used to compute the permanent can 
also be used to compute the determinant, by appropriately selecting the signs. 


We recall now the dehnition of the parity function, and we extend it to ordered 
partitions. 


Definition 3.1. Let D, Y be ordered sets of the same size. For a bijection vr : D —)■ 
Y we dehne its sign or parity as sgn( 7 r) := (—where N{'k) is its number of 
inversions: 


N{n) := |{(a,a') E : a < a', 7r(a) > 7r(a')}| . 

Let y = (Yi,..., Yk) be an ordered partition of Y, i.e., Y = Yi U ■ ■ ■ U Y),. We dehne its 
sign to be sgn(3^) := where N{y) is: 

^(3^) := lUVi^yj) ■ Vi e Yi, yj E Yj, i < j, yi > yj}\. 

Equivalently, we can associate to 3^ a permutation vr^ : {1, 2,..., |y |} —)■ Y that consists 
of blocks: we put hrst Yi (sorted), then Y 2 (sorted), and so on. Then sgn(3^) = sgn(7r^). 

From the dehnition above it is clear that we can obtain the sign of a permutation 
in O(n^) by counting the number of inversions. However, it is well known that we can 
hnd it in 0{n) by counting its cycles. 

For a matrix M, there is a natural order for its row set A and column set X, namely 
from top to bottom and from left to right. We recall the dehnition of the determinant: 

Det(M) = ^sgn(7r) ][[[ Ma,^(a) 

aeA 


TT 
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where the sum is over all bijections tt : A ^ X. Similarly, for a fixed matrix M and for 
some D A and Y O X we define the partial determinants: 

(6) det{D, Y) := ^ sgn(7r) JJ Ma,n{a) 

TT aG-D 

where the sum is over all bijections tt : D ^Y. Note that Det(M) = det{A,X). 

We now provide a recursion formula similar to the one in Lemma We need one 
lemma before. 


Lemma 6. Let D, Y he ordered sets, and let tt : D Y be a bijection, which we view 
as a subset of D x Y . Let V = {Di, ..., D^) and y = (Yi,..., Y^) be partitions of D 
and Y . Let vr = tti U ■ ■ ■ U he a decomposition with Tij C Dj x Yj. Then 

k 

sgn(7r) = sgn('D)sgn(37) JJsgn(7rj). 

i=i 


Proof. It follows from the multiplicativity of the sign function. 


□ 


Lemma 7. Under the same conditions of Lemma\^ then 


(7) det{AT,,Y) = sgn{V)^sgn{y)det{At,Yt) JJ det{AT,.,Yc^ 

y i=i 


where det{-, ■) is as in 0 the sum is over all y = (Y*, Yci, ■ ■ •, b'c*) satisfying Q, 
V = {At,AT ,^, 


A 


and sgn(-) is as in Definition 


3.1 


Proof. The proof is basically the same as the one of Lemma The only difference is 
that we have the additional factor sgn(7r), but it factors because of Lemma □ 


Despite the resemblance between equations (|^ and ([^, the latter is not a subset con¬ 
volution because of the sign factors. Therefore, we cannot use the algorithm from 
in this case. We now show to the complexity analysis. 


Lemma 8. Given the values of the partial determinants det{At,Yt) and det{AT^,,YcA, 
we can evaluate equation Q for all Y satisfying ([^ in 0{k{n + 3‘^)). 

Proof. We will hrst express equation Q in a similar format as formula ([^ of Lemma 
To simplify the notation, let y =: (Yq, Yi, ..., Y^). For each j let Yf = Yq U Yi U ■ ■ ■ UTj-, 
and observe that sgn(3^) = sgn(YQ^~^, Tj). Then equation ([^ can be rewritten as: 

k 

D(Y) = sgn(P) Y. n L) D,Y 

You-uYfc=yi=o 

where Yq, ..., W C y(t) and 

A := x{Tt) \ X{t) Ao := 0 A, := \ x{t) 

D{Y) := det{AT,, Y U A), Do(Yo) := det{At, Yo U Aq), D,{Yj) := det{AT^.,Yj U A,) 

Yff = YoU---U Yj A^ := Ao U ■ ■ ■ U Aj 
S,{Yf-\ Y,) := sgn(Yo'-' U Ar\ Y,- U A,) 
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For each 0 < / < A;, and for each Yq C x(A), let 

i 

= sgn(®) ^ ]l^i{yt'.yi)D,,{yi) 

YbU-UYi=Y(* j=0 

Note that Dq(Y) = D(Y), and thus it is enough to compute for all 1. We can do 
this recursively, observing that = sgn('D)Zi)o(?) and 

(8) £>J+'(FJ+') = 5^ s,+i(y„‘,y.+i)cJ(yJ)A+i(r,+,) 

Y({uy+i=Yj+i 

We reduced the problem to evaluating the above formula, and we will show that for 
each I we can do this in 0{n + 3^^). Assume for now that the signs S';+i(yJ, FJ+i) are 
known. Then for each of cardinality i, we can evaluate ([^ in 0(2*). Thus, for 
all Yq'^^ we require - (“)2*) = 0(3"'), where w = \x{t)\- We will see that after 
a precomputation that takes 0(n), we can obtain 5i+i(io,^+i) 1^ which will 

complete the proof. 

Observe that 

S'z+i(Fo', Yi+i) = sgn(Fo, ^z+i) sgn(F((, Az+i) sgn(A|), Yi+i) sgn(A(), Az+i). 

Note that the last factor does not depend on the partition and it can be precomputed 
in 0(n). Also note that the hrst factor can be computed in 0{u) = 0(1), so we can 
ignore it. We are left with the second and third factor. 

For each x G x(t), let 

= \{y e Az+i : x>y}\. 

We can precompute for all x in 0{ujn) = 0{n). Note that sgn(yj, A^+i) = (—1)'^ 

where N = YlxeY • Thus, after the precomputation, we can obtain this factor in 

0{u) = 0(1). A similar procedure can be done for sgn(AQ, Fz+i). □ 

Theorem 9. Let M be a matrix with associated column graph . Let {T,x) be a tree 
decomposition of of width uj. Then we can compute Det(M) in 0{n^ + n4"'). 

Proof. There are two changes with respect to the proof of Theorem]^ First, in the base 
case we need to compute the determinant of all submatrices of Mq. Using an expansion 
by minors as in the proof of Lemma we can do this in 0(4"*), i.e., the same as for the 
permanent. Second, for the recursion formula we use LemmaThis increases the time 
per node from 0{kt 2‘^) to 0{kt{n+ 3"^)). Therefore, the overall cost is 0{n‘^ + n4P). □ 

We conclude this section by presenting an open question. Given the resemblance 
in the dehnition of the permanent and the determinant it is not surprising that they 
can be computed using very similar tree decomposition methods. The immanant of a 
matrix is another closely related notion: 

ImmA(M) := E Xa('^) Ma^7r{a) 

7T a£A 
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where the sum is over all bijections tt ■. A ^ X and xa is an irreducible character 
of the symmetric group. The immanant reduces to the permanent when is the 
trivial character, and it reduces to the determinant when xa is the sign character. The 
computational complexity of immanants has been analyzed in e.g., 13 . A natural 


question that arises is whether a tree decomposition method can be used to compute 
them. We remark that the recursion in Q does not hold for the immanant as xa is not 
necessarily multiplicative. 


Question. Given a matrix M of bounded treewidth, can we compute ImmA(M) in 
polynomial time? In particular, can this be done if either the height or the width of the 
Young diagram is bounded? 


4. Bipartite decompositions 


In the previous section we showed a decomposition method based on the column 
graph . We showed that we can compute the permanent in 0{nY^^), where ux is 
the treewidth of G^. In this section we will extend this decomposition method to work 
with the bipartite graph G (see Dehnition 2.2). We will show that we can compute the 
permanent in 0{n2^), where u is the treewidth of G. A Matlab implementation of our 
algorithm is available in www.niit.edu/~diegcif, 

Let G be the bipartite graph of M. As in the previous sections, we index the rows 
with a set A and the columns with X. We now rephrase the dehnition of a tree 
decomposition of G. A bipartite decomposition of G is a tuple (T, Q!,x), where T is a 
rooted tree and a : T —)■ {0,1}^, x : ^ {0,1}^ assign some a(t) C A and x(^) ^ Ai 

to each node t of T, that satishes the following conditions. 


i-1. The union of is the whole row set A. 

i-2. The union of {x(t)}ier is the whole column set X. 

ii. For every edge (a, x) of G there exists a node t of T with a E a(t),x G x(t)- 
iii-1. For every a E A the set {t : a E a(t)} forms a subtree of T. 

iii-2. For every x E X the set {t : x E x{t)} forms a subtree of T. 

The width oj of the decomposition is the largest of \a{t) \ + |x(^)| among all nodes t. 
Note that the above literals are consistent with the ones in Dehnition 12.11 

As before, we now present an example to illustrate the use of the bipartite graph for 
computing the permanent. 



Figure 4. Bipartite graph G of the matrix of Figure]^ and a tree de¬ 
composition T. 
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Example 4.1. Consider again the matrix M of Figure]^ Note that Perm(M) can also 
be expressed in the following form: 

Perm(M) = perm({ai, 02}, {xi, X4}) perm({a3, 04, 05}, {x2, X3, X5}) 

+ perm({ai, 02, 03}, {xi, X3, X4}) perm{{ai, 05}, {x2, X5}) 

- Ma3,x3 perm({ai, 02}, {xi, X4}) perm({a4, 05}, {x2, X5}) 
perm({ai, 02, 03}, {xi, X3, X4}) = perm({ai, 02}, {xi, X4}) + Ma^^x^ perm({oi, 02}, {xi, X3}) 

perm({o3, 04, 05}, {x2, X3, X5}) = Ma^^x^ perm{{aA, 05}, {x2, X5}) + perm{{ai, 05}, {X3, X5}) 

To evaluate the above formula we need to compute four 2x2 permanents, and we need 
in total 16 multiplications. It turns out that this formula arises by considering the tree 
decomposition of the bipartite graph shown in Figure]^ 

4.1. Bipartite decomposition algorithm. Algorithm presents our dynamic pro¬ 
gramming method to compute Perm(M) using a bipartite decomposition. As for Al¬ 
gorithm [TJ for each node t we compute a table Pt, following a topological ordering of 
the tree. The permanent of M is in the table corresponding to the root. There are two 
main routines: SubPerm computes the permanents of all submatrices, and EvalRe- 
CURSION evaluates a recursion formula, which is slightly more complex than the one of 
Algorithm [1} 

As before, the values Pt{D,Y) computed correspond to a partial permanent of the 
matrix. Consider the collection 

5 = {{D, Y) : a(TO \ a{t) CDC a(T,), x(T,) \ x(t) C F C y(T,)}. 

Observe that {D, F) G S' is completely determined by {D fl a{t),Y fl x(t)). Therefore, 
if we let .D = D n a(f), F := F O x(f), there is a one to one correspondence between S, 
and the collection S := {{D,Y) : D C a{t),Y C x(f)}. The partial permanents that 
we are interested in are 

Pt{D, Y) := perm{D, Y) = perm{D U ia{Tt) \ a{t)), Y U (xiTt) \ x{t)))- 
For the root node Proot{c>i{root), x{root)) = perm{A, X) = Perm(M). 

4.2. Recursion formula. The recursion formula that method EvalRecursion of 
Algorithm evaluates is given in the following lemma. 


Lemma 10. Let M be a matrix with associated bipartite graph G. Let {T,a,x) be a 
bipartite decomposition of G. Let t be an internal node of T , let Tt Y T denote the 
subtree rooted in t, and let D, Y be such that 

(9) a{Tt) \ ait) CDC a(Ti), xiTt) \ xit) C F C y(Tt), \D\ = |F| 

Let tc^,... the children oft. Then 

k 

(10) permiD,Y) = y^perm{Dt,Yt) ]^(-l)I^F lperm(Ac,-, ^cc,) 

v,y j=i 
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Algorithm 2 Permanent with bipartite decomposition 


Input: Matrix M and tree decomposition {T,a,x) of bipartite graph G{M) 
Output: Permanent of M 
1: procedure BipartPerm(M, T, a, x) 

2 : order := topological ordering of T starting from its leaves 

3 : for t in order do 

4: Qt := SUBPERMS(f, M) 

5 : if f is a leaf then 

6: Pt := Qt 

7: else 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


Cl,..., Cfc := children of t 
Pt :=EVALRECURSION(f, Qt, Pci, • • • , PcQ 
return Proot{c({root), xiroot)) 
procedure SuBPERMs(f, M) 

Qt{D, Y) ■= perm{D, Y) for all D C a{t),Y C x(t) 

procedure EvALRECURSlON(f, Qt, Pc^, ■ ■ ■, Pc^) 
for Cj child of t do 

A" := a{cj) \ a(t). A" := a{cj) fl a(t) 

■= X{cj) \ xiti AJ := xicj) n x{t) 

QtcjiD,?) := {-l)\^\Qt{D,Y) for all D C A“,y C Aj 
Qcc, (P, Y) := Pc, (P U A“, Y U Aj) for all D C A“, Y C Aj 

Pt := SUBSETCONVOLUTION(Qi, Qtci,Qcc^, • ■ ■ , Qtc^^QccJ 
procedure SubsetConvolution(Po, Pi, , P2k) 

P(P, Y) := Po(Po, Fo) Pi(Pi, ?,)■■■ P2k{D2k, Ffc) 


DoU---UD2k=P 

YoU-UY2k=Y 


where perm{-, ■) is as in Q and the sum is over all V = {Dt, Dtc^,...), y = (Yt, Ytc^, ■ ■ ■) 
satisfying: 


D = DtU (Pici U Pcci U ■ ■ ■ U Dtc, U PccJ 
Y = YtU (Eici U Ecci U • • • U Ytc, U Y^Q 
a{Tc.) \a{t) C Pec, P a(Pci) Dt C a{t) Dtc^ C a(t) n a(fc,) 

a(^c,) \ xit) P Ycc^ Y x(Tc,) Yt C x{t) Ytc^ C x{t) H xiPj)- 

To prove this lemma we need some additional notation. We view a bijection vr : P — 
F as a subset of P x F, by identifying it with the set {(a,7r(a)) : a G P}. For a given 
node t and for some P, F satisfying (|^, we denote 

perm*{D,Y) : = EH 

TT a<^D 






19 


where the sum is over all bijections n : D ^ Y such that 

(11) TT n (a(t) X x(t)) = 0- 

We now show a different recursion formula, which is closer to the one in Lemma 


Lemma 11. Following the same notation as above, the following eguation holds: 

k 

(12) perm{D,Y) = ^ perm{Dt,Yt)Y\_perm*{Dc.,Yc.) 

Yt,Ycj,Dt,D,^ 3=1 

where the sum is over all Yt,Yc., Df, such that 
(13a) D = DtU (Hci U ■ ■ ■ U 

(13b) F = FtU(y;, u-.-uwj 

(13c) a(^c,) \ a(t) C C a{Tc.) Dt C a{t) 

(13d) x(Te,) \ X{t) C W, C x{Tc,) D C x{t). 

Proof. Let tt : H —>■ F be a matching, which we view as a subset oi D xY. Note that 

D = {Dn a{t)) U (L) n a(TcJ \ a{t)) U ■ ■ ■ U {D n a{TcJ \ a{t)) 

Y = {Yn x{t)) u (F n x(TeJ \ x{t)) u ■ ■ ■ U (F n x(TeJ \ x{t)) 

Let’s decompose tt in a similar way as above. Let nt be the submatching of tt with 
domain contained in H fl a{t) and range contained in F fl xit)- Equivalently, vr* = 
71 n {D D a{t)) X (F n x(t)). Observe that if some a G H ft Q!(t) is not in the domain of 
7it, then a G a(Tc), 7 r(a) G x(Ec) for some child c. The reason is that by dehnition of tree 
decomposition, there must be some node tp with 7 r(a) G xifp)^^, G o.{tp). However, the 
assumption on a says that 7 r(a) G x(Tt)\x(t) and thus we must have tp G Tc for some c. 
Similarly, if some a; G F fl x(t) is not in the range of iit, then x G x(Tc), G a(Tc) 

for some child c. In the same way, we have the following 

7i{D n a{Tc) \ a{t)) C DO xiTf) 

7r-i(Fnx(T,)\x(t))CFna(T,) 

The above paragraph implies that we can decompose 

(14) TT = TTt U TTci U • ■ ■ U TTcj, 


in such a way that Tit F (D H a(f)) x (F fl x(t)) and for each child c we have that 
TTc ^ (-D n a{Tc)) X (F n x(Tc)) and tTc satishes (0. Moreover, it is easy to see that 
such decomposition is unique. 

Let Fj C F n xit) be the range of Tit and let Fc C F fl x(Ec) be the range of tTc. 
Analogously, dehne Dt F Dna{t) and Dc F ZlnQ!(Tc) as the domains of Tit, t^c- Observe 
that x(Ec) \ x(t) ^ Ec, as if a: G x(Tc) \ xif) then (7r“^(a:), a:) G Tic by construction, and 
thus X G Fc. Similarly, x(Tc) \ x(t) C Fc. Therefore, the equations (13) are satisfied. 

Thus, for any matching ti Y D x Y there is a unique partition as in (14) such that 
the ranges and domains Yt,Yc, Dt, Dc satisfy (13). On the other hand, assume that 
Yt,Yc, Dt, Dc satisfy (13), and we are given some matchings 7it,7ic with these ranges 
and domains and such that 0 holds. Then equations ([I^ tell us that we can merge 
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them into a matching vr with domain D and range Y. Condition (0 ensures we are 
not overcounting, as it implies that decomposition (14) is unique. These remarks imply 
equation ( 12 ). □ 


Remark 4.1. Let {T^,x) be a tree decomposition of the column graph , and let 
{T,a,x) be the corresponding decomposition of the bipartite graph G given in Exam¬ 
ple 2^ In such case, the above lemma reduces to Lemma 


We now derive the recursion formula in Lemma [10 
lemma by using inclusion-exclusion. 


which follows from the above 


Proof of Lemma 10 


For a child c of t, we will show that 


(15) perm*{Dc,Yf) = ^ {-iy^^‘^^perm{Dtc,Yi,f)perm{Dcc,Ycc)- 

DtcUDcc=Dt, 

YtcUYcc=Yc 

DtcCa(t),YtcCxit) 


Combining equations (12) and (15) we obtain equation (10), concluding the proof. 
Given a matching tTc : He —t Y^, let /(tTc) := tTc fl {a{t) x x(f)) and let Ia{T^c) Y a{t), 
Y x(t) be the domain and range of /(tTc). For some Dtc Y DcHait), Ytc C YeP^xif) 
with I Ac I = IA I, let 


perm*{DcYx, Dtc, Ytc) := 

ttc-Dc^Yc a 
ioL {j^c)—Y)tc 

Note that perm* {Dc,Yc) = perm*{Dc,Yc',^,^). Observe now that given matchings : 
Dtc —t Ytc and tTcc : Dc\Dtc —t Yc\Ytc, we can merge them into a matching 7i* : Dc ^ Yc 
that satisfies — Dtc, -^x(^c) — Therefore, we have the following equation 

perm{Dtc, Ytc) perm{Dc \ D^, W \ Y^) = E perm*{Dc, W; Dtc Yc)- 

Dtc^Dtc 

yt*D.ytc 


Based on the above formula, we can now find perm* { Dc , Yc ) using inclusion-exclusion 
(or Mobius inversion): 


perm*{Dc,Yc-,%,%) = E E {—If perm{Dtc, Ytc) perm{Dc \ Dtc, Yc \ Ytc). 

i DtcQDcna{t) 

YtcCYcCixit) 

\DtcHYtc\=i 


Rewriting the above equation leads to (15), as wanted. 


□ 


4.3. Complexity analysis. We just derived the recursion formula (10) which is used 
in Algorithm As in the proof of Lemma [^ this formula is a subset convolution and 
thus it can be evaluated efficiently using the algorithm from 
time of Algorithm is stated now. 


10 . The overall running 


Theorem 12. Let M be a matrix with associated bipartite graph G. Let {T,a,x) be a 
bipartite decomposition of G of width u. Then we can compute Perm(M) in 0(n2‘^). 
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Proof. The proof is very similar to the one of Theorem For each node f G T, we 
will compute perm{D^Y) for every pair D,Y that satishes equation ([^. We will show 
that for each t we can compute perm{D,Y) for all such D,Y in 0{{kt + 1)2'^), where 
kt is the number of children of t. Observe that for the root node tr we will compute 
Perm(M) = perm{a{Tr), x{Tr))■ This will conclude the proof. 

The base case is when f is a leaf of T, so that Tt = {t}. Let Mq be the submatrix of M 
with rows a{t) and columns x(t). We need to obtain the permanent of all submatrices 
of Mq. As |a(f)| + |x(^)| < we can do this in 0(2‘^) using Lemma [i 

Assume now that t is an internal node of T with kt children and let T*, T" that 
satisfy (9). Then equation (10) tells us how to hnd perm{D,Y). Similarly as in 
Lemma 4, we can evaluate this formula in 0{kt2‘^), assuming we know the values 
of the terms in the recursion. Note that we already found perm{Dcc-, Ycc) for all chil¬ 
dren in the recursion. We can hnd perm{Dt,Yt) for all Dt,Yt in 0(2‘^) in the same 
way as for the base case, and this includes the values perm(Dtc,Ytc). Then, it takes 
0(2^ -|- kt 2^) = 0{{kt + 1)2^^) to compute perm{D, Y) for a all D, Y. □ 


Similarly as in Theorem we can hnd an analogous algorithm for the determinant. 

Theorem 13. Let M he a matrix with assoeiated bipartite graph G. Let {T,a,x) be a 
bipartite decomposition of G of width u. Then we can compute Det(M) in Oipn? +niY). 


Proof. We just need to follow the steps of Section 3.5 For instance, the recursion is 


det{D,Y) = Y,^gn{V)sg-n{y)det{Dt,Yt) det{Dtc^,Yt,^)det{D,,^,Y,,^) 

v,y j=i 


where y = {Yt,Ytc,Ycc),V = {Dt,Dtc, Dec). The complexity analysis is basically the 
same as in the proof of Theorem The base case can be done in 0{2‘^) using an 
expansion by minors. The recursion can be evaluated in 0{kt{n + 3‘^)) in a similar way 
as in Lemma [HI □ 


5. Mixed discriminant and higher dimensions 

The mixed discriminant of n matrices is a common generalization of the permanent 
and the determinant. As such, it is also hard to compute in the general case. We show 
now that the techniques presented earlier generalize to compute mixed discriminants. 
Even more, we show that this method extends to compute similar functions in higher 
dimensional tensors. 

5.1. Mixed discriminant. Let M be a list of n matrices of size n x n. Equivalently, 
we can think ofMasanxnxn array. We index the hrst coordinate with a set A, 
and the second and third coordinates with sets X^,X‘^. The mixed discriminant of M 
is given by 

Disc(M) := ^ sgn(7r^)sgn(7r^) JJ 

TT^.TT^ a^A 
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where the sum is over all bijections tt^ : A —)• and tt^ : A —)■ and sgn is 

the parity function. For a G A, let Ma denote the n x n matrix obtained by fixing 
the hrst coordinate. Observe that if Ma = m for some matrix m and for all a E A, 
then Disc(M) = n!Det(m). In the case that Ma is diagonal for all a E A, then 
Disc(M) = Perm(Zi)) where D is the matrix obtained by concatenating these diagonals. 
Some of the properties of mixed discriminant are discussed in [^. 

In the case of a n x n matrix, a bipartite graph was the natural structure to represent 
its sparsity. Similarly, if we are given a sparse n x n x n array M, a natural structure 
is a tripartite graph G, as follows. Let G be the graph on A U U X^, where for each 
nonzero entry Ma^x^_,x 2 P^t a triangle {a, a;i,X 2 }. 

We rephrase the dehnition of a tree decomposition of a tripartite graph G. A tripartite 
decomposition of G is a tuple (T, a, x^, where T is a rooted tree, a : T —)■ {0,1}"^, 
: T —{0,1}^ and '■ T —>■ {0,1}^ , that satisfies the following conditions. 

i. The union of {a{t)}t£T (resp. x^ X^) is the whole A (resp. X^, X^). 

ii. For every triangle {a,xi,X 2 ) in G there is a t with {a,xi,X 2 ) G (a x x x^)(f)- 

iii. For every a E A (resp. X^, X^) the set {t ■. a E a{t)} is a subtree of T. 

The width of the decomposition is the largest of \a(t)\ + |x^(f)| + \x^(t)\ among all 
nodes t. Note that the above literals are consistent with the ones in Dehnition 12.11 In 
particular, observe that the second condition does not impose additional constraints 
due to Lemma [U 

We proceed to extend the previous results to the mixed discriminant. For some sets 
D CA,Y^ GX^ and C X^ we denote 

(16) disc{D, Y\ Y^) := ^ sgn(7r^)sgn(7r^) JJ M^i^a)nHa) 

TT^jTT^ aGD 

where the sum is over all bijections 7i^ : D ^ Y^ and vr^ : D —>■ This only makes 

sense if \D\ = |X^| = |X^|, and otherwise we can dehne disc{D,Y^,Y‘^) = 0. 

As for the case of the permanent, the dynamic program to compute Disc(M) has two 
main steps: computing the mixed discriminant of all subarrays of M, and evaluating 
some recursion formula. For the hrst step, it is easy to see that the approach from 
Lemma extends, as we show now. 

Lemma 14. Let Mq be a ni x n 2 x array. Let Ao,Xq,Xq he its set of coordinates, 
and let S = {{D,Y\Y'^) C Aq x X^ x X^ ; |D| = IX^I = \Y^\}. We can compute 
disc{D, Y^, Y^) for all triples in S in 0{n^ax where Umax = max{ni, 77 , 2 , n^}- 

Proof. For i = 1, 2,..., min{ni, n 2 , n^} let 

Si := {{D,Y\Y^) C Ao X Xl x X" : \D\ = \Y^\ = \Y‘^\ = t}. 

We can hnd disc{D, X^, X^) for all {D, X^, X^) G Si using the values of Si_i as follows. 
Let oq be the hrst element in D, it is easy to see that 

disc{D,Y^,Y‘^) = ^ e{xi,X 2 )Mao,xi,x 2 disc{D\ao,Y^ \xi,Y^ \x 2 ) 

Xl&Y^,X2e:Y'^ 
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where e{xi,X 2 ) is either +1 or —1. To be concrete, if we identify with the set 

,i}, then e(ji, j 2 ) = (—Thus, for each triple (D, Y^, y^) we just need to 
loop over 7121 x 3 terms, and for each we need O^rimax) to hnd D \ oq, \ xi, \ X 2 . □ 

The recursion formula we need to evaluate is given in the following lemma. 

Lemma 15. Let M he a list of n matrices of size nxn, with associated tripartite graph 
G. Let {T,a,x^,X^) be a tripartite decomposition ofG. Let t be an internal node ofT, 
let Tt YT denote the subtree rooted in t, and let D, F^, F^ be such that 


a{Tt) \ a{t) Y D C a{Tt), 


\D\ = IF^I = |F^ 


(17a) 

(17b) x\T,) \ x\t) YY^Y x\T,), x\T,) \ x\t) YY^ Y x\T,) 

Let Cl,... ,Ck be the children of t. Then 

- 2 ^ 


(18) 


disc{D,Y^,Y^) = ^ sgn{y^)sgn{y^)disc{Dt,Yf,Yf 


n(-i)l°“J F4, y4)d«c(c„„ ry, y; 

i=i 


V = (D„ A„. ..= {yfxx ...), = (vy Yx ... 


Z = Z, U (Z,„ u z„, u ■ ■ ■ u z„. u z, 


where disci-, is as in (16), sgn(-, •) as in Definition 3.1, and the sum is over all 

^ ^,2 _ ^v'2 ^.2 Satisfying: 

where Z e {D,Y\Y^} 

((Tcj) \C{t)Y Zcc^ Y ({Tc.) Zt Y ({{) Ztcj Y({t)n ((4,) where ( G {a, x\ X^}- 

Proof. Let vr^ : D —F^ and vr^ : H —)■ F^ be matchings. Observe that Lemma says 
that sgn(7r^) will factor in the tree decomposition, leading to the term sgn('D)sgn(3^^). 
Similarly, sgn(7r^) leads to the term sgn(Il)sgn(3^^) (note that sgn(Il) cancels). There¬ 
fore, for the rest of the proof we can ignore all sign factors. We can think of the pair 


(vr^, TT^) as a subset of H x F^ x F^. In a similar way as we did in the proof of Lemma 11 
there is a unique decomposition of (vr^, tt^) of the form 

(19) (ttS vt^) = (ttI, TTt) u (7r^\, TT^J U ■ ■ ■ U (vr^\, TT^J 

where (7r^^,7r^) C (a x x y^)(t), and for each child c we have that (vr^,7r^) F (a x 
X x^){Tc) and 

(20) {nl. Til) n (a X y^ X x^W) = 0- 

Note that by construction tt^ , vr^ have the same domain. Let Dt denote this domain, 
and let Yf,Y^ denote the respective ranges. Similarly, let Dc,Yf,Y^ be the domain 
and ranges of vr^, vr^. Then we have 

(21a) Z = ZtU (Zci U ■ ■ ■ U Fcj.) where Z G {D, Y^, F^} 

(21b) ((Tc) \Cit) YZcY ((Tc) Zt Y ((t) where C e {a, y\ y^} 


Thus, for matchings (tt^, tt^) G 71 x F^ x F^ there is a unique partition as in (19) and 
the corresponding domains and ranges satisfy (21). On the other hand, assume that 
Dt, Dc, Yf,Yf, Y^, Yf' satisfy (13), and we are given some matchings 7 ij, 7 il, 7 i^, vr^ with 
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these domains and ranges and such that (vr^,vrc) satisfy (20). Then equations (21) tell 
us that we can merge them into matchings 7r^,7r^ with domain D and ranges 


Condition (20) ensures we are not overcounting. Then we have 


( 22 ) 


dtsc{D,Y\Y^) = 




J=1 


where the sum is over all triples as in (21), and where disc*{Dc-,yc-iyc) is similar to 
disc{Dc,Y^,Y^), except that it only uses matchings ( 71 ^ 5 '^c) satisfying (20). 


Finally, we can obtain disc*{Dc, Y^, Y^) using inclusion-exclusion in a similar way as 


in the proof of Lemma IT 
(23) 




f-l)IO“l Y,i Y,l)disc(D„ ri, Yl). 


DtcUDcc=Dc, DtcQait) 

YguY,l=Y,^, YgCxHt) 


Combining equations (22) and (23), we obtain equation (18). 


□ 


We proceed to the complexity analysis. 


Theorem 16. Let M he a list of matrices with associated tripartite graph G. Let 
(T, he a tripartite decomposition of G of width u. Then we can compute 

Disc(M) in 0{n^ Y nS"^). 


Proof. The proof is very similar to the one of Theorem 12 For each node f G T, we 


compute disc{D, Y^, Y^) for every triple D, Y^, Y"^ that satisfies equation (17). We will 
show that for each t we can get disc{D, Y^, Y'^) for all such D, Y‘^ in 0{{kt -|- l)(n-|- 
3^)), where kt is the number of children of t. 

The base case is when t is a leaf of T. Let Mq be the subarray of M given by indices 
{a(t), X^(t)). Then we need to find the mixed discriminant of all subarrays of Mq. 

We can do this in 0(2'^) using of Lemma 


14 


Assume now that t is an internal node of T with kf children and let D,Y^,Y^ that 
satisfy (^. Then equation (18) tells us how to hnd disc{D,Y^,Y‘^). Similarly as in 
Lemma q we can evaluate this formula in 0{kt{n Y 3“^)), assuming we know all terms 
in the recursion. We already found disc{Dcc, Yf,, Y^) for all children in the recursion, 
and we can find disc{Dt, Yf, Yf') for all Dt, Yf, Y^ in 0(2'^) in the same way as for the 
base case. This leads to a bound of 0{{kt -I- l)(n -|- 3^^)) to compute disc{D, Y^, for 
all D,Y\Y‘^. □ 


5.2. Higher dimensions. It is easy to see that our methods extend to compute gen¬ 
eralizations of the permanent and determinant in higher dimensions. We consider a 
square (d-l- l)-dimensional array (or tensor) M of length n, i.e., of size nx ■ ■ ■ xn {dYl 
times). Here we assume d to be constant. Let’s index the first coordinate of M with 
a set A, and the following coordinates with sets X ^,..., X'^. Consider a function F of 
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the form 
(24) 


F{M) = E n ei(7ri) ■ ■ ■ €d{7ld)Ma,^i(a),...,iTdia) 

Q'GA 


where the sum is over all bijections tti : A ^ X\ and where ei{'n'i) is either 1 or sgn(7r;). 

Let’s consider two special cases of equation ( [24| ). The simplest case is when ei{Tri) = 1 
for all 1. We refer to such function as the {d+ l)-dimensional permanent and we denote 
it as Perm(M) [^. Some applications of this permanent are shown in 

Consider now the case when d+1 is even and ei{7ii) = sgn(7r/) for all /. This is perhaps 
the simplest generalization of the determinant, and it is usually referred to as the first 
Cayley hyperdeterminant 14 . Some applications of the hyperdeterminant are shown 
in 


As opposed to the 2-dimensional case, computing the hyperdeterminant is 


25 . 

As 

hard 



We now proceed to extend our decomposition methods to this setting. We asso¬ 
ciate a. {d + l)-partite graph G where for each nonzero entry of M we put a {d + 1)- 
clique in the respective coordinates. A tree decomposition of G can be seen as a tuple 


{T,a,x\ 


• • ,X 


The width oj of the decomposition is the largest of |Q!(f)| -f- |x^(t)| + 


-|- x“(t) among all nodes t. 


As before, for some sets D ^ A,Y’' O X\ we consider the function 
(25) f{D,Y\. 


E ) :— 'y ^ ■ ■ ■ ed{T^d)Ma^-ni{a),...,-Ka{a) 

7ri,...,7rcf a£A 


where the sum is over all bijections ni : D ^ YK 

There are two steps in order to generalize our results to this setting: evaluate / in 
all subarrays, and evaluate the recursion formula. For the former, the approach from 
Lemma 1^ (and Lemma 14) has a simple generalization. Indeed, Barvinok shows this for 
the case of the hyperdeterminant [^. The proof is the same for an arbitrary function 
F as in (|24|). Thus, we have the following. 


Proposition 17 ( J^)- Let Mq he a {d Y 1)- dimensional array of size no x • • • x Ud, 
and let f be as in ( [^ . Let Aq, Xq, ..., Xq be its set of coordinates, and let 

s := {(T),y\...,y'^) c Ao X X X : \D\ = |y^| = •■■ = 

We can compute f{D, Y ^,..., Y'^) for all tuples in S in 0(n((^)^ where Umax = 

max{no, ...,nd]. 


Repeating the same analysis as in the proof of Lemma lA, the recursion formula is: 
f{D,Y, S,{V,y')---UV,y‘)f(D„Yp..,Y,'‘) 

k 

j](-i)i°“ji/(A=,, i'4,..., .., 14 ) 

i=i 


(26) 
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where Si{V,y^) is either 1 or sgn('D)sgn(V) depending on ei, and the sum is over all 
tuples V,y^,... ,y^ such that 

Z = ZtU [Ztc-, u Zcci U ■ ■ ■ U Ztc^ U ZccJ where Z G {D, 

C{Tcj) \ C{t) Y C ({T^.) Zt c ({t) Ztc^ C ((t) n C(tcJ where C e {«, • • •, 

The complexity of the decomposition algorithm is as follows. 


Theorem 18. Let M he a square {d + 1)-dimensional array of length n, with {d + 
l)-partite graph G. Let F be a generalized determinant/permanent as in (24). Let 
(T, a, x^, ■ ■ ■, x‘^) be a tree decomposition of G of width oj. Then we can compute F{M) 
in 0{nf + n 3^). 

Proof. The proof is very similar to past ones. For each node t, we compute f{D, ..., 
for all valid tuples. We show that for each t we can do this in 0{{kt + l)(n + 3‘^)), where 
kt is the number of children of t. 

The base case, i.e., leaf nodes, reduces to Proposition!^ leading to a bound of 0(2^). 
For an internal node t, equation (26) tells us how to find f{D, ..., Y’^). Similarly 
as in Lemmawe can evaluate this formula in 0{kt{n + 3‘^)). □ 

For the special case of the permanent, we can give a better bound. 


Theorem 19. Let M be a square [d + 1)-dimensional array of length n, with (d + 1)- 
partite graph G. Let (T, a, ..., x‘^) be a tree decomposition of G of width oj. Then 
we can compute Perm(M) in 0(n2‘^). 

Proof. If there are no sign factors we can follow the procedure of Lemma for the 
recursion, leading to a bound of 0{{kt + 1) 2‘^) per node. □ 


6. Mixed volume of zonotopes 


The mixed volume MVol of n convex bodies Ki, ..., Kn in M” is the unique real 
function that satishes the following properties. 

• MVol is multilinear and symmetric in its arguments. 

• MVol(iF,..., K) = n\ vol{K), where vol denotes the volume. 

Alternatively, it can be shown that the function /(A) := vol{'^.XiKi) for Aj > 0, is a 
homogeneous polynomial, and MVol is the coefficient of Ai ■ ■ ■ A„. For more information 
about mixed volumes, see e.g., 30 . We focus here in the case that all bodies Ki are 


zonotopes, which are a special class of polytopes. 


6.1. Mixed volumes and permanents. 

Definition 6.1. A zonotope is a polytope that is a Minkowski sum of line segments, 
i.e., it has the form 

z = [0, l]zi + [0, 1 ]Z 2 H-h [0, l\zm = {riZi H-h : 0 < r* < 1} 

where Zi G M” are vectors. In case zi, ... ,Zm are linearly independent, we say that is 
a parallelotope. 

The mixed volume of zonotopes has a simple description as follows. 
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Proposition 20. Let ® zonotope, for i = 1,... ,n. Then 

(27) 


Jn ^ 


MVol(z\ ..., z’") = ^ |Det(4,4’ • • • ’ 

Jl V ■ ijn £ '^n 

Proof. The multilinearity of the mixed volume implies 

MVol(z\ ..., z") = MVol([0,1]4,..., [0,1]4J. 

Thus, we just need to argue that 

MVoldO, 1]4,..., |0,1|4J = |Det(4,..., 4)1 

which follows by noting that ^ parallelepiped with sides XiZf, and thus 

its volume is given by the absolute volume of the determinant. □ 


The mixed volume of n parallelotopes reduces to a permanent when their main axes 
are aligned, as shown now. 


Corollary 21. Let ui,... ,Un G ML and let M G ]R>q” ® nonnegative matrix. Let 

z'^ = ® zonotope. Then 

MVo1(2:^, ..., z^) = |Det(Mi,..., Un) \ Perm(M). 


Proof. We just need to use equation (27), and cancel out all terms that contain a re¬ 
peated vector Uj, as the determinant is zero. The remaining terms have |Det(Mi,..., | 

as a factor and we get the desired formula. □ 


6.2. Graph representation. To use a decomposition method for mixed volumes we 
need to have a graph description of the zonotopes. We consider now two different 
graphs that can be associated to a set of zonotopes, and more generally to polytopes. 
The hrst one is a bipartite graph that can be thought of as the analogue of the bipartite 
graph of a matrix. The second one has to do with the sparsity in the standard basis 
representation and it is a more intuitive notion for general polytopes. 

Definition 6.2. Let Q be a set of n polytopes in M^. Let U denote the set of all vectors 
(up to scaling) that are parallel to some edge in Q. We refer to U as the edge directions 
of Q. The edge graph G{Q) is a bipartite graph with vertices QUU and edges (g,u) if 
q contains an edge parallel to u. 

Definition 6.3. Let Q be a set of n polytopes in M^. Let X = {xi,..., denote the 
coordinates. The coordinates graph G^{Q) has X as vertex set, and for each polytope 
q E Q we form a clique in all its non constant components. Note that if g = 
is a zonotope, we form a clique in the nonzero coordinates of zj. 

Remark 6.1. Note that the edge graph G{Q) is invariant under affine transformations 
of Q, whereas the coordinates graph G^{Q) is not. 
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Example 6.1 (Zonotopes of bounded treewidth). Consider the following zonotopes: 
(28a) = [0, l](aie„ + ei) + [0, l]ei 

(28b) z* = [0,1] (ojCn + Cj - Ci-i) + [0,1] (cj - ei-i) for i = 2,..., n - 1 

(28c) Z'^ = [0, ijianCn - Cn-l) 

where {ej}j is the canonical basis and ai E Tj are some integers. 

Note that the segments of the above zonotopes are all nonparallel (there are 2n + 1 
edge directions). This means that the edge graph G has n connected components, one 
for each of the zonotopes, and each of component is either a 2-path or a 1-path. Thus, 
G has treewidth 1. As for the coordinates graph G^, it is the union of the triangles: 
Xi := {xj, Xj+i, a;„}. It is easy to see that G^ has treewidth 2. 


The following example shows the relationship between these graphs and the matrix 
graphs we used before. 

Example 6.2 (Relationship with matrix graphs). Let tie zonotopes 


as in Corollary 21 The edge graph G of the zonotopes has vertices Z VJ U where 


Z = {z'^}i and U = {uj}j, and it has an edge {z\uj) whenever Mij ^ 0. If we replace 
Z with the row set and IJ with the column set, this is precisely the bipartite graph of 
matrix M. 

On the other hand, the coordinates graph G^ depends on the sparsity structure of 
vectors U. Assume now that U = {ej}j is the canonical basis. Then G^ has an edge 
{xj, Xk) whenever there is some z^ with Mjj ^ 0 and Mj ^ ^ 0. This corresponds to the 
column graph of M. 


Because of the above example it is expected that the tree decomposition methods 
derived for the permanent should allow to compute mixed volumes of certain families 
of zonotopes. Indeed, we show now if there are few edge directions and the edge graph 
G has bounded treewidth then the mixed volume can be computed efficiently. 


Theorem 22. Let Z he a set ofn zonotopes in W^. Let U be the set of edge directions of 
Z, and assume that d := \U\ — n is constant. Let G denote the edge graph of Z. Given 
a tree decomposition of G of width u, we can compute MVol(Z) in 2‘^). 


Proof. If \U\ < n it follows from (27) that MVol(Z) = 0. If \U\ = n, then Corollary 21 
tells us that we just need to compute the determinant of the ufs and the permanent 
of some matrix M. We can find the determinant in 0{n^) with linear algebra. For the 
permanent, as the edge graph G corresponds to the bipartite graph of M, we can use 
Theorem |l2[ Thus, we can hnd this permanent in 0{n2‘^). 

If \U\ > n, let W U of size n. For each z E Z assume 2 ; = Xluef/lO’ 
some scalars c„,. Let zw = l]cnM, and let Zw = {zw '■ z E Z}. It follows from 


equation (27) that 


MVol(Z) = 


WCU, \W\=n 


MVo\{Zw). 


For each such W the associated graph is a subgraph of G, so we can compute MVol(Z^) 
in 0{n^ -|- n2‘^). As there are = 0{n'^) possible W, the result follows. □ 
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Example 6.3 (Zonotopes with few edge directions). Consider the following zonotopes: 


z* = [0, IjaiCi + [0, l]bie 


for i = 1,... ,n 


where {ej}j is the canonical basis, e := Xlj = (I 5 • • •) 1) and ai,bi G Z are some 
integers. Observe that there are n + 1 edge directions: ei,..., e„, e. Also note the that 
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the edge graph G is a tree, consisting of pairs ( 2 ;*, e*) and ( 2 ;*, e). Therefore, Theorem 
says that can compute the mixed volume of these polytopes in polynomial time. 

As an application, recall that Bernstein’s Theorem gives a correspondence between 
mixed volumes and the roots of systems of polynomials. This allows us to conclude that 
we can efficiently count the number of solutions of the following system of equations: 


0 = c. 


:,i + Ci,2 + ^,3 + 0^,4 2:“' 


X, 


for i = 1,..., n. 


6.3. Hardness result. Theorem 22 shows that it is possible to exploit tree decompo¬ 
sitions for mixed volume computations. However, it restricts the zonotopes to have a 
small number of edge directions. This is a strong requirement which is not satisfied 


in many cases (such as in Example 6.1). Unfortunately, we will see that we need this 


condition. We remark that the same condition appears in discrete optimization, where 
it allows to derive (strongly) polynomial time algorithms for certain discrete convex 
optimization problems [^ . 

We show now that computing the mixed volume of the zonotopes in Example 6T is 
^(^P-hard. This shows that mixed volumes of zonotopes continue to be hard, even if 
both G and have bounded treewidth. We use a similar reduction as in 
they prove that the volume of zonotopes is ^P-hard. 


19 , where 


Lemma 23. The determinant of the following n x n matrix is si -1- S 2 + • • • + 


M = 


r 1 -1 0 0 ■■■ 0 0-1 

0 1 -1 0 ■■■ 0 0 

0 0 1 -1 ■■■ 0 0 


0 0 0 0 ■■■ -10 

0 0 0 0 ■■■ 1-1 

-.Sl *2 .S 3 *4 ■■■ Sn - 


Proof. If we perform Gaussian elimination we end up with an upper triangular matrix 
where the diagonal is: = 1 for i < n and = si s„. □ 


Proposition 24. The following problem is ffP-hard. Given integers oi,..., On, com¬ 
pute the mixed volume of the n zonotopes of eguations (28). 


Proof. We consider the ^^^P-complete problem Subset-Sum: given a set of integers A, 
determine the number of subsets S' C A with sum zero. Let A: = n — 1, A = {oi,..., Ok) 
and let On = S he a. parameter. We will show that the solution to the Subset-Sum 
problem is given by 

(29) ^MVol(z\ ..., zt-i) - MVol(z\ ..., zto) + ^MVol(z\ ..., z^i)- 
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Let’s evaluate equation (27). Consider the following n x (2n — 1) matrix. 


1 

1 

-1 

-1 

0 

0 ■ 

0 

0 

0 

0 

0 1 

0 

0 

1 

1 

-1 

-1 ■ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 ■ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 ■ 

■ -1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 ■ 

1 

1 

-1 

-1 

0 

0 

0 

0 

0 

0 

0 ■ 

0 

0 

1 

1 

-1 

_ ai 

0 

d2 

0 

«3 

0 ■ 

■ a*;-l 

0 


0 

5 


Ms = 


Observe that columns 2i — 1 and 2i of correspond to z*, and the last column 


corresponds to z^. Then formula (27) considers submatrices of that use columns 
ji ,... ,jn where jt G {2i — 1, 2i} for i = 1,..., A; and = 2n — 1. Note now that 
for any subset S' C A, there is a natural submatrix to consider: if a* G S then 

ji = 2i — l and otherwise ji = 2i. This correspondence is a bijection. Observe also that 
each submatrix has the form of Lemma Thus, we have the following equation: 




SCA 


6 ^ ^ di 


Finally, observe that for any integer s we have 


1 .N , , , 1 m , 1 if s = 0 

2 '^ ^ ' ' ' 2' ' lO otherwise 


The last two equations imply that (29) indeed counts the subsets S with sum zero. □ 
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